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numerical,  and  experimental  investigations  are  presented.  Introducing  a  constriction 
into  an  annular  resonator  splits  each  longitudinal  duct  mode  into  two  modes,  one  of  a 
higher  fi'equency  with  a  pressure  antinode  at  the  constriction  and  one  at  a  lower 
fi’equency  with  a  velocity  antinode  near  the  constriction.  The  lower  mode  is  more 
sensitive  to  changes  in  the  length  and  porosity  of  the  constriction  than  the  higher  mode. 
Overall  agreement  between  measured  and  predicted  mode  shapes  and  resonance 
fi-equencies  is  veiy  good.  It  was  found  that  it  is  necessary  to  include  end 
corrections  at  the  constriction  to  get  accurate  agreement  between  measured  and 
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INTRODUCTION 


The  purpose  of  this  thesis  is  to  investigate  acoustic  resonances  in  a  constricted, 
annular  resonator.  The  motivation  for  this  project  is  based  on  a  related  problem  in  the 
area  of  thermoacoustis  -  the  annular  prime  mover.  The  intent  is  that  the  measurements 
reported  in  this  thesis  will  serve  as  preliminary  validation  of  models  being  developed  for 
the  prime  mover  project. 

The  a:Hmiitha1  modes  of  a  uniform  cross  section  annulus  are  degenerate.  The 
introduction  of  a  constriction  splits  the  degeneracy  resulting  in  two  modes  having  different 
resonance  frequencies  and  distributions  of  acoustic  pressure  (mode  shape).  The  goal  of 
this  thesis  is  to  both  measure  and  predict  the  frequency  splitting  and  mode  shapes. 

A  theoretical  discussion  will  be  covered  in  the  Chapter  H.  Chapter  HI  is  a 
description  of  the  experimental  apparatus  including  the  annular  acoustic  resonator  and 
peripheral  equipment.  In  addition,  a  detailed  discussion  of  the  data  collection  procedure 
will  be  covered.  Chapter  IV  is  an  analysis  of  the  data  obtained  and  comparison  with 
predictions.  The  thesis  is  summarized  and  important  conclusions  are  addressed  in  Chapter 
V. 


1 


2 


n.  THEORY 


A.  INTRODUCTION 

The  purpose  of  this  chapter  is  to  outline  the  theoretical  analysis  of  a  driven, 
constricted  annular  resonator.  The  discussion  begins  with  the  solution  of  the  wave 
equation  in  and  the  effective  radius  of  a  uniform  cross  section  annular  resonator.  This 
solution  is  developed  in  cylindrical  coordinates.  However,  because  the  purpose  of 
this  thesis  is  to  investigate  frequencies  below  those  required  to  excite  cross  modes,  we 
adopt  a  simplified  plane  model  for  the  majority  of  analysis. 

The  next  step  is  to  introduce  a  constriction.  It  is  shown  that  the  constriction 
results  in  frequency  splitting  of  what  were  degenerated,  orthogonal  modes  in  the  uniform 
cross  section  resonator.  It  is  also  shown  that  the  degree  to  which  a  particular  mode  is 
excited  depends  upon  the  location  of  the  driver  relative  to  the  constriction.  Because  of 
dissipation,  the  modes  overlap  in  frequency  space  and,  in  general,  two  modes  are 
simultaneously  excited.  When  this  situation  occurs  the  resonance  frequencies  of  the 
individual  modes  cannot  be  determined  from  peaks  in  the  frequency  response. 
However,  pole-zeros  analysis  of  the  frequency  response  can  yield  accurate  values  for  the 
resonance  frequencies.  The  chapter  concludes  with  a  discussion  of  pole-zero  analysis  of 
two  coupled  modes. 

B.  SOLUTION  FOR  A  UNIFORM  ANNULAR  RESONATOR 


The  analysis  begins  with  the  general  solution  of  the  a(X>ustic  wave  equation 
in  cylindrical  coordinates 


where 


V 


1  d 


r  dr 


^  d 

r  - 


t) 


dt  ’ 

1  d  ' 
r  "  de 


d  " 
dz^ 


(2.1) 

(2.2) 
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Assuming  harmonic  time  dependence,  the  wave  equation  becomes  the  Helmholtz 
equation 


V^P  = -k^V 


where  = 


-Y 

<cJ 


(2.3) 


The  general  solution  of  Eq.  (2.3)  is  assumed  to  be  expressed  as  the  product  of  three 
terms,  each  a  function  of  only  one  variable. 


Figure  2-1.  Geometry  for  the  annular  resonator. 


?(r,0,z)  =  R(r)&(0)Ziz) 

Substituting  Eq.  (2.4)  into  (2.3)  and  dividing  by  RQZ  ,  gives 


1  1  d 

(  dR  " 

\  \  \  d^®  \  d^Z  2 

R  r  dr  \  dr  > 

Rearranging  to  isolate  the  z 

\  \  d  (  dR\ 

)  ^  r  ^  0  de""  '  Z  dz'^ 

-  dependent  terms  gives 

1  1  ^'0  .  ,  \  d^z 

R  r  dr  ^ 

J  dr) 

r"  0  Z  dz^ 

Following  the  standard  separation  of  variables  solution  technique , 


(2.4) 


(2.5) 


(2.6) 


let 


1 


=  -  k 


Z  dz'^ 

The  general  solution  to  this  equation  is 

Z  =  Aj  cos(A:jr)+  sin(^jz) 

Equation  (2.6)  now  becomes  after  multiplying  through  by 


(2.7) 


(2.8) 

and 


rearrangmg 


4 


_r _ £_ 

R  dr 


Let 


1  d'^% 


ir  =  - 


m 


1  d^% 
0  dd^ 


©  dO^ 

The  general  solution  to  this  equation  is  conveniently  expressed  in  the  form 

©  (^  )  =  CO  s(m  0  ) 

Equation  (2.9)  is  now  reduced  to 


r  d 
R  dr 


[r  '-/«")=  0  , 


where 


-  k] 


(2.9) 

(2.10) 

(2.11) 

(2.12) 

(2.13) 


The  general  solution  to  this  equation  is 


(2.14) 


where  jJ^r^  and  Y^{k^7^  are  the  cylindrical  Bessel  and  Neumann  functions.  The 
Neumann  function  must  be  retained  because  the  origin  is  excluded  in  our  annular 
resonator. 

All  boundaries  are  assumed  to  be  rigid,  so  the  normal  component  of  acoustic 
velocity  must  be  zero  at  the  boundaries.  Applying  rigid  boundary  conditions  at  z  =  0 
and  z  =  implies  that  =  0  and 

;=0,1,2...  (2.15) 

The  function  0(0)  must  be  periodical,  in  other  words,  0(0)  =  0(0  +  2;r) .  This 
condition  is  satisfied  if  m  is  an  integer,  m  =  0,1,2  ••••. 

Application  of  rigid  boundary  condition  at  r  =  aand  r  =  b  implies  that 
h„J'„{k^a)  +  B„Y\„{k^a)  =  0 
A„J'».(M)  +  B„r'».(M)=  0.  (2.16) 

Expressed  in  matrix  form  these  equations  become 
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'j'Sko) 

U.M  icMJvbJ  UJ-  ^  -  > 

The  values  of  k,a  that  satisfy  this  equation  are  those  for  which  the  determinant  of  the 


first  matrix  is  zero ,  i.e. 


J'„(k,a)Y‘„(k^b)  -  J  „.{k,b)Y'„{k,a)  =  0  .  (2.18) 

Let  the  values  of  k^a  that  satisfy  Eq.  (2.18)  be  called  a„„ .  With  this  definition,  the 


first  of  Eqs.  (2.16)  becomes 


=  -A, 


Y'ACC„n)  ’ 


(2.19) 


and  Eq.  (2. 14)  becomes 


^(r)= 


(2.20) 


For  the  purposes  of  this  thesis,  we  are  interested  in  fi'equencies  low  compared 
to  those  required  to  excite  cross  modes.  To  eliminate  pressure  variations  in  the  z 
direction,  we  require  that  /  =  0  in  Eq.  (2.15).  In  this  case  Z(z)  =  1,  and  k^  =  k^ . 
We  are  also  interested  in  modes  for  which  the  circumference  of  the  resonator  is 
approximately  one  wavelength  long,  i.e.  m  =  1  in  Eq.  (2.11).  Under  these  conditions, 
the  solution  to  the  Helmholtz  equation  is 


P(r,0)  =  A  {J,{k,r)  -  -^4^^Yi(V)]cos(0  +  <!>). 

\^\n) 


(2.21) 


The  dimensions  for  our  particular  resonator  are  a—  10.00  cm  and  b  —  15.30  cm.  The 
lowest  value  of  is  then  =  0.796  .  Recall  that  for  /=0, 


Also,  k  =  — ^ .  Combining  these  expressions  for  k^ 
a 


implies  that 


2;r 

X  a 


(2.22) 


Let  us  define  the  effective  circumference  of  the  resonator  —  m  X  ,  ox  for  m  —  1 

=  X  .  Combining  this  definition  with  the  Eq.  (2.22)  and  rearranging  yields 
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2  na 


(2.23) 


Alternatively,  we  can  define  an  effective  radius  ^  ^  ■  Using  Eq.  (2.23)  ,  we 

have  .  For  =  0.796  ,  R^ff  =  12.56  cm. 
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C.  PLANE  WAVE  THEORY 


The  parameter  R^ff  is  defined  in  such  a  way  that 


R  =  hjH--  ^  — 

In  In  k  ' 

The  angular  coordinate  0  can  be  related  to  the  arc  length  x  measured  at 


(2.24) 
R^ff  through 


e  =  =  kx  .  (2.25) 

With  this  definition ,  the  0  -  dependence  of  the  acoustic  pressure  becomes 

cos(^x  +  ^).  (2.26) 

This  function  can  be  expressed  in  terms  of  complex  exponentials  as 

cos(kx  +  ^)  =  (2.27) 

Substituting  this  result  into  Eq.  (2.21)  gives 

P(r,^,0=  .  (2.28) 

The  interpretation  of  Eq.  (2.28)  is  that  the  pressure  is  represented  by  two  counter 
propagating  waves  that  behave  essentially  like  plane  waves.  This  analogy  will  be 
used  to  analyze  the  nature  of  the  sound  field  in  the  annular  resonator. 


The  plane-wave-like  solution  will  be  valid  for  fi-equencies  below  those  required 
to  excite  the  first  cross  mode.  The  frequencies  corresponding  to  normal  modes  in 
the  r-zplane  are  given  by 
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(2.29) 


f.  =  +  * 


27t  a 


Table  2-1  lists  a  few  values  of  a, 


m\n 

0 

1 

2 

1 

0.796 

6.025 

11.903 

2 

1.589 

6.119 

11.987 

3 

2.378 

6.465 

12.124 

Table  2-1.  A  few  values  of  a  . 


In  our  annulus  a  =  0.10  m,  so  it  is  seen  that  is  approximately  60  for  m  = 


7t  (X 

1.  Similarly,  4  is  0.05  m,  so  —  is  equal  to  62.8,  or  nearly  equal  to 

Lt. 


Therefore  the  frequencies  of  the  first  cross  mode  in  either  the  z  or  r  direction  are 
approximately  equal  and  have  value  in  the  vicinity  of  3.3  -  3.6  kHz.  Our  primary 

c 

interest  is  in  frequencies  corresponding  to  =  A.,  or  f  =— —  =  440  Hz. 


1.  Uniform  Cross  Section  Annular  Resonator 
a.  Undriven  case 

As  explained  above,  for  low  frequencies  the  armular  resonator  can  be 
modeled  as  a  straight  pipe  of  length  with  periodic  boundary  conditions  at  the 

two  end  as  shown  in  Fig.  2-2.  Assume  that  there  are  two  traveling  waves  in  pipe. 
One  is  propagating  in  the  positive  x-direction  and  represented  by 

.  (2.30) 

The  other  propagates  in  the  negative  x- direction  and  given  by 

P_  =  (2.31) 
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The  total  pressure  in  the  pipe  is  a  sum  of  these  two  pressures  and  given  by 

P,  =  +  P_  =  .  (2.32) 

These  subscripts  get  used  later  in  a  different  way  use  p  and  n  for  positive  and  negative 
instead  of  1  and  2. 


Figure  2-2.  Annular  resonator  and  equivalent  strm^t  resonator. 


Also,  the  volume  velocity  U  can  be  obtained  from  knowledge  of  the  pressure  and  the 


o  c 

pipe’s  acoustic  impedance,  Z  =  — where  S  is  cross  section  area  of  pipe, 

j3 

Sk 


U,=U^-U  =  — . 


(op, 

Periodic  boundary  conditions  are 

P  =  P 
and 

where  L  =  L  .  Substituting  Eq.  (2.32)  into  (2.34)  gives 

A  =  B- 


(2.33) 

(2.34) 

(2.35) 

(2.36) 


1-e-'^  ■ 

Substituting  this  expression  into  Eq.  (2.33)  and  applying  Eq.  (2.35)  yields 

cosAZ.  =  1 .  (2.37) 


This  implies  that 


kL  =  2m7i 


(m  =  0,l,2 . ), 


(2.38) 
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It  should  be  pointed  out  that  has  slightly  different  values  for  different  values  of 
m. 

b.  Driven  case 

In  addition  to  knowing  the  resonance  frequencies  of  the  annular 
resonator,  we  are  also  interested  in  the  mode  shapes,  i.e.  the  acoustic  pressure  and 
velocity  as  function  of  angular  position.  There  is  no  unique  orientation  of  the  modes 
in  an  undriven,  uniform  cross  section  annulus.  Either  the  system  must  be  driven,  or 
a  nonuniformity  in  cross  section  must  be  introduced  to  fix  the  mode  shapes.  In  the 
undriven  uniform  cross  section  case,  any  superposition  of  two  degenerate 
orthogonal  modes  (  one  proportional  to  cosmO  ,  the  other  to  sinm0  )  satisfies  the 
boundary  conditions.  In  order  to  find  the  mode  shape,  it  is  necessary  to  install  a 
driver  in  resonator.  Assume  that  a  branch  driver  is  located  at  x  =  0  as  shown  on  Fig 
2-3.  The  volume  velocity  injected  by  the  driver  is  given  by 

(2.40) 

x=0  x=L 


Figure  2-3.  Driver  installed  straight  resonator. 

The  boundary  condition  for  pressure  is  given  by 

P  =P 

which  is  same  as  that  for  the  undriven  system.  This  conditions  leads  to 

A(l-c-'“)  +  B(l-e'“)  =  0. 


(2.41) 

(2.42) 


The  boundary  condition  for  volume  velocity  now  includes  an  additional  term  and  is  given 
by 

(2.43) 


which  yields  the  equation 


A(l-e-'“)-B(l  +  e*)  =  ^«„. 
Equations.  (2.42)  and  (2.44)  can  be  expressed  in  matrix  form  as 


'l  - 

'A' 

0 

_1  - 

-  1  + 

B 

Po^ 

^  U  X) 

(2.44) 


(2.45) 


This  equation  has  the  form  MC  =  U .  The  value  of  the  coefficients  A ,  B  can  be 
determined  from  C  =  M~^U .  Once  A  and  B  are  known,  the  pressure  and  volume 
velocity  are  found  from  Eqs.  (2.32)  and  (2.33). 


2,  Constricted  Annular  Resonator. 

a.  Frequency  vs  area  ratio  (no  driver) 

The  analysis  of  the  constricted  resonator  is  more  complicated  due  to  its 
non  uniform  cross  sectional  area.  The  geometry  is  shown  in  Fig.  2-4. 


x=0  x=  / 


x=L 


Figure  2-4.  Constricted  resonator  with  no  driver. 


Calculation  of  the  resonance  frequency  of  a  constricted  resonator  is  similar  to  that 
of  the  uniform  cross  section  resonator.  In  order  to  solve  this  case,  the  basic  equations 
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of  pressure  and  volume  velocity  in  each  section  have  same  form  as  the  uniform  cross 
section  resonator. 


P,  =  , 

(2.46) 

(2.47) 

At  the  junction  of  two  sections,  pressure  and  volume  velocity  are 

continuous  which  gives 

the  following  equations 

P,(x  =  /)  =  P2(x  =  /), 

(2.48) 

P,(x  =  0)  =  P2(x  =  Z), 

(2.49) 

U,{x  =  t)  =  U,(x  =  D, 

(2.50) 

U,{x  =  0)  =  U^ix  =  L). 

(2.51) 

These  equations  can  be  combined  into  a  4  x  4  matrix  in  the  following  form. 

gjV  _  g-J'M  _  1 

1  _  ^ 

.  ®  ^0.  (2.52) 

^  5,Jt,  S,k^  C 

_1  -  S';/:;  ^-Jk,L  ^2^2  ^Jk,L  D 

This  matrix  has  the  form  of  MC  =  0  .  The  resonance  frequencies  are  those 
frequencies  that  satisfy  the  equation  det[M]  =  0. 

Unlike  the  uniform  cross  section  case,  the  introduction  of  a  constriction 
splits  the  degeneracy  of  the  modes.  The  two  modes,  called  the  high  frequency  and 
low  frequency  modes,  correspond  to  different  mode  shapes.  ( In  what  follows,  the 
word  “frequency”  be  will  frequency  dropped  and  the  modes  referred  to,  simply,  as 
“high”  and  “low”.)  The  high  mode  has  a  pressure  antinode  centered  on  the  constriction 
whereas  the  low  mode  has  a  velocity  antinode  centered  on  the  constriction.  The 
difference  in  frequency  between  the  high  and  low  mode  increases  with  increasing 
length  or  decreasing  cross  sectional  area  of  the  constriction. 
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The  dependence  of  the  frequency  splitting  is  shown  in  Figs.  2.5  a-c. 
The  plots  show  resonance  frequency  as  a  function  of  area  ratio  for  three  different 
constriction  lengths.  The  resonator  contains  air  at  20®C  and  atmosphere  pressure.  Area 
ratio  is  defined  as  the  ratio  of  the  cross  sectional  (open)  area  in  the  constriction  to 
that  of  the  resonator.  In  the  limits  of  zero  area  ratio,  i.e.  when  the  constriction  is 
completely  blocked,  the  high  mode  frequency  corresponds  to  A  = 

and  the  low  mode  frequency  corresponds  to  X  =  2(L^  ~  ^ctmanction)  •  These  limits  are 
appropriate  for  the  two  lowest  modes  of  a  rigid,  rigid  resonator  of  length  ^constriction  * 


calculated  resonance  freq  vs  area  ratio 


Figure  2-5a.  Resonance  frequencies  with  a  long  constriction 


16  J 
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b.  Mode  shapes  of  a  constricted  resonator 

Calculation  of  the  pressure  distribution  in  a  driven,  constricted  resonator 
is  also  similar  to  that  for  the  driven  uniform  cross  section  case.  For  the  driven,  constricted 
resonator,  the  resonator  is  divided  into  three  sections,  as  shown  in  Fig.  2-6.  The  driver 
is  located  at  x  =  d. 


x=d 

Figure  2-6.  Constricted  resonator  with  driver. 


This  system  gives  rise  to  the  following  set  of  equations  for  pressure  and  volume  velocity  in 
the  three  regions: 


P,  = 

(2.53) 

(2,54) 

P3  = 

(2.55) 

’  op,  (op. 

(2.56) 

^2^2  ^^jXait-k.yx')  ^2^2  £j^ /((W+M) 

"  cop,  (op. 

(2.57) 

^  ^2^2  ^^i(<ot-k,x)  ^2^2  p^K<ot*k^x) 

^  cop,  cop. 

(2.58) 

(2.59) 

The  three  different  sections  also  give  rise  to  six  boundary  conditions: 
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P,(jc  =  0)  =  P3(x  =  L), 

(2.60) 

U,{x  =  0)  =  U,{x  =  L), 

(2.61) 

P,(x  =  /)  =  P3(x  =  /), 

(2.62) 

ti 

il 

11 

(2.63) 

(2.64) 

=  </)  +  Mo  =  . 

(2.65) 

Application  of  these  boundaiy  conditions  to  Eqs.  (2.53)  -  (2.59)  yields  a  set  of  equations 


wMch  can  be  combined  into  a  matrix  equation  of  the  form  MC  =  U . 


1 

1 

0 

0 

! 

_ 1 

'A' 

B 

0 

1 

-1 

0 

0 

-^2^  -ik^L 

SA, 

^2^2  jk^L 

SA, 

0 

r\ 

0 

0 

C 

0 

<^2^2  ^Jk^l 

0 

0 

D 

0 

SA 

SA 

p 

0 

0 

0 

J-j 

0 

0 

-e^ 

F_ 

^  Ud 

(2.66) 


Solving  the  inverse  form  of  Eq.  (2.66)  gives  the  value  of  the  coefficients  A  through  F. 
Substitution  of  these  coefficients  into  Eq.  (2.53)-(2.55)  gives  the  mode  shapes. 

(1)  Mode  Selection  and  Driver  Location.  Because  of 
dissipation,  the  quality  factor  (Q)  of  a  particular  mode  is  finite.  This  means  that,  in 
general,  when  the  resonator  is  driven  at  a  fi'equency  in  the  vicinity  of  the  resonance 
of  either  the  high  or  low  mode,  both  modes  will  be  excited.  In  order  to  measure 
the  mode  shape  of  a  single  mode,  it  is  necessary  to  excite  only  that  mode  as  nearly  as 
possible.  This  can  be  accomplished  to  an  acceptable  degree  by  proper  location  of 
the  driver.  If  the  driver  is  positioned  at  a  pressure  antinode  of  a  particular  mode 
and  driven  at  the  resonance  fi'equency  of  that  mode,  then  that  mode  is  the 
predominant  one  excited. 

Mode  selection  is  illustrated  the  next  sequence  of  figures  which  show  the 
calculated  frequency  response  and  mode  shapes  for  the  resonator  with  a  constriction 
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length  of  and  area  ratio  0.7.  For  Fig.  2-7,  the  driver  is  located  180®  from  the 

center  of  the  constriction  which  should  correspond  to  a  pressure  antinode  for  the 
Wgh  mode.  Figure  2-7a  shows  the  frequenqr  response  as  indicated  by  computed 
magnitude  of  acoustic  pressure  measured  45®  from  constriction.  This  figure  shows  that 
there  is  only  a  single  peak  in  the  frequency  response.  The  peak  frequency,  439.0  Hz, 

corresponds  to  the  resonance  frequency  of  the  high  mode  for  a  long  constriction 

with  a  0.7  area  ratio,  as  shown  in  Fig.  2-5b.  Figure  2-7b  shows  the  calculated  mode 
shape  when  the  resonator  is  driven  at  439.0  Hz.  It  is  seen  that  there  is  indeed  a  pressure 
anitnode  at  both  the  constriction  and  driver  location.  It  is  also  evident  that  the  standing 
wave  ratio  is  very  high,  another  indication  of  single-mode  excitation. 


Frequency  Response 


Figure  2-7a.  Frequency  response  of  the  constricted  resonator  when  driver  is  located  180® 
from  the  constriction  and  driven  at  the  high  mode  resonance  frequency  for  the 


length,  0.7  area  ratio  constriction. 
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Mode  Shape  for  Length  L/32,  Area  Ratio  0.7,  and  f=439.0 


Figure  2-7b.  Mode  shape  when  the  resonator  is  driven  at  the  peak  frequency  as 
determined  from  Fig.  2-7a. 

Figure  2-8a  shows  the  frequency  response  when  the  driver  is  moved 
to  a  location  90®  from  the  constriction.  This  driver  location  should  be  near  a 
pressure  node  for  the  high  mode  and  pressure  antinode  for  the  low  mode.  The 
frequency  response  again  shows  only  one  peak.  Comparison  with  Fig.  2-5b  indicates 
that  this  mode  is  the  low  mode.  Figure  2-8b  shows  the  mode  shape.  It  is  seen  that 
now  there  is  a  pressure  node  at  the  constriction.  Again,  the  standing  wave  ratio  is 
high  indicating  single-mode  excitation. 
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Frequency  Response 


Figure  2-8a.  Frequency  response  of  the  constricted  resonator  when  driver  is  located  90° 
from  the  constriction  and  driven  at  the  low  mode  resonance  frequency  for  the 
length,  0.7  area  ratio  constriction. 


Mode  Shape  for  Length  L/32,  Area  Ratio  0.7,  and  f=427.2 


Figure  2-8b.  Mode  shape  when  the  resonator  is  driven  at  the  peak  frequency  as 
determined  from  Fig.  2-8a. 
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In  the  cases  of  the  driver  located  at  180  degrees  or  90  degrees  from 
the  constriction,  the  driver  excites  only  one  mode,  either  high  or  low.  If  the  driver  is 
placed  at  a  position  other  than  180  degree  or  90  degree,  two  modes  will  be  excited 
simultaneously.  Figure  2-9  shows  the  frequency  response  when  the  driver  is  located  45 

degrees  from  the  constriction.  We  have  used  the  0.1  area  ratio  constriction  to 

make  the  amount  of  frequency  splitting  larger. 


Frequency  Response 


ratio  constriction. 

We  can  see  from  Fig.  2-9a  that  both  the  high  and  low  modes  are  excited.  In  this  case, 
they  are  well  separated,  in  other  words,  there  is  no  significant  overlap.  When  there 
is  little  overlap,  it  is  easy  to  pick  up  the  resonance  frequency  of  each  mode.  A  constriction 

of  small  area  ratio  (0.1)  and  long  size  22.5  degree),  as  in  Fig.  2-9a,  separates  the 

two  modes  well  enough  for  there  to  be  little  overlap.  Figure  2.9b  shows  the  mode 
shape  when  the  resonator  is  driver  at  the  lower  peak  frequency.  This  figure  shows 
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that  even  when  the  driver  is  at  a  location  which  would  allow  both  modes  to  be 
excited,  when  the  modes  are  well  separated,  single  mode  excitation  is  still  possible. 


Mode  Shape  for  Length  L/16,  Area  Ratio  0.1,  and  f=320.0 


Figure  2-9b.  Mode  shape  when  the  resonator  is  driven  at  the  lower  of  the  two 
peak  frequencies  from  Fig.  2-9a. 


Figure  2- 10a  shows  the  frequency/'  response  for  the 


length,  0.7  area  ratio  constriction  used  in  earlier  figures.  In  contrast  to  Fig.  2-9a,  the 
modes  show  significant  overlap.  This  overlap  causes  two  problems.  First,  it  means 
that  a  driver  located  at  45°  will  excite  both  modes  to  a  significant  degree,  even 
when  driver  at  one  of  the  peak  frequencies.  Second,  the  peak  frequencies  no  longer 
correspond  to  the  resonance  frequencies.  Overlapping  is  most  significant  for  large  area 
ratio,  short  constrictions. 

Figure  2- 10b  shows  the  mode  shape  when  the  resonator  is  driven 
at  the  lower  peak  frequency  in  Fig.  2- 10a.  It  is  seen  that  the  standing  ratio  is  low  in 
this  case,  indicating  excitation  of  more  than  one  mode.  The  next  section  addresses  a 
technique  to  determine  resonance  frequencies  in  the  case  of  significant  overlap. 
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Figure  2- 10a.  Overlapped  wave  form  of  frequency  response  with  driver  45  degrees  from  a 
constriction  which  has  a  0.7  area  ratio  and  >  (5.625  degree)  length. 


Mode  Shape  for  Length  L/32,  Area  Ratio  0.7,  and  f=428.5 


Figure  2- 10b.  Mode  shape  when  the  resonator  is  excited  at  the  lower  peak  frequency 


D. 


POLE -ZERO  ANALYSIS 


For  understanding  the  acoustics  of  an  annular  resonator,  it  is  important  to  find 
accurate  values  of  the  resonance  fi’equencies.  These  resonance  fi’equencies  are 
independent  of  the  driver  location  but  are  dependent  upon  what  kind  of  constriction  is 
being  used  -  the  length  and  area  ratio.  As  we  have  seen,  the  driver  may  wtcite  one  or 
two  modes  simultaneously  depending  upon  its  excitation  fi’equency  and  location.  When 
only  one  mode  is  exdted,  it  is  simple  to  find  resonance  fi'equency  by  finding  the 
fi'equency  corresponding  to  the  peak  in  the  fi'equency  response  .  But  for  the  case  when 
the  driver  excites  two  modes  simultaneously,  the  peak  fi’equency  cannot  be  used  to 
determine  the  resonance  fi’equency.  This  situation  was  discussed  in  the  previous  section. 
Pole-zero  analysis  can  be  used  to  determine  the  resonance  fi’equencies  in  the  case  of 
significant  modal  overlap. 


1.  Analysis 


For  a  driven  simple  harmonic  oscillator,  the  complex  speed  is  represented  by  [Ref 


1] 


Fe ' 


U  = 


R.  +  -  £l 

The  expression  for  u  can  be  written  as 


(2.67) 


U 


jG)  — e 
m 


jatt 


(2.68) 


oj^nij 


by  introducing  the  resonance  fi'equency  and  quality  factor,  Q  =  ° /r  ■ 

By  analogy,  we  assume  the  pressure  in  a  resonator  driven  near  resonance,  at 
one  point  in  space,  has  the  form 
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(2.69) 


where  Q  is  now  the  quality  fector  for  the  resonator.  Factoring  the  denominator  of  Eq. 
(2.69)  results  in 


This  equation  can  be  simplified  if  the  denominator  is  written  as  [Ref  2] 

(©  +  jP\ X®  +  ./>2 ) .  (2.71) 


where  the  poles  and  are, 


(2.72) 


where  Q>  i  =  to  „  ,  for  the  single  mode  case.  The  poles  are  written  in  the  form  {a -\-jp) 
because  the  Hewlett  Packard  signal  analyzer  used  for  experimental  measurement  outputs 
the  poles  in  that  form.  Equation  (2.74)  can  be  simplified  into  a  form  given  by 
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For  two  overlapping  modes,  the  complex  pressure  amplitude  is  the  sum  of  the  individual 
amplitudes 


P.,  =  P,  +  P: 


(2.76) 


and  has  the  form 

p  _  Ai(m  +yp4)  + AaCm  H-y/?;) 

(<y+y>i  ) . (®  +  7>4) 


(2.77) 


So,  even  though  the  two  modes  are  overlapped  and  the  amplitude  of  the  individual 
responses  are  added,  the  poles  still  show  separable  frequencies. 

The  numerator  of  Eq.  (2.77)  can  be  factored  in  the  same  way  as  the  denominator 
of  Eq.  (2.69).  In  this  case  the  numerator  is 


where 


a  + 


-  jc  +  d 


2(Aj  +  A^)) 


0)  + 


^  -Jc-d  ^ 


2(Ai  +  Aj). 


(2.78) 


c  =  A,(-|-).A.(-|;-), 


(2.79) 


and  d  = 


A, 


02 


+  A, 


0. 


+  4(a,  +  A,  Xa,®,’  +  Aj®,=  ) .  (2.80) 


From  Eq.  (2.78)  we  can  find  zeros. 

The  resonance  frequency  and  quality  factor  can  be  obtained  from  the  poles  using 
Eq.  (2.75).  The  resonance  frequency  is  then 

G)  ^  =  -s/a  ^  (2.81) 
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and  the  quality  factor  is 


2a  ' 


(2.82) 


If  Q  is  large,  an  approximation  can  be  used  to  simplify  Eq.  (2.81).  The  term  of 


can  be  neglected  since 


r  ^  * 


K2QJ 


is  small  compared  to  o  .  So,  b  = 


For  data  analysis,  we  fit  the  measured  fi-equency  response  using  pole-zero 
analysis.  The  fit  parameters  give  the  resonance  frequency  and  Q  of  the  individual  modes. 


E.  THERMO-VISCOUS  DAMPING 


Up  to  this  point,  losses  have  been  ignored.  However  any  real  gas  has  viscosity 
and  will  exchange  heat  with  the  walls  of  the  resonator.  The  boundary  condition  at 
the  walls  ofthe  resonator  is  that  the  normal  velocity  is  zero.  When  viscosity  is  included 
the  correct  boundary  condition  has  both  the  normal  and  tangential  velocity  going 
to  zero.  This  creates  a  viscous  boundary  layer  which  contributes  wscous  losses  to  the 
propagation  of  sound  [Ref.  3,4].  Similarly,  the  gas  near  the  walls  will  exchange 
heat  with  the  walls.  This  thermal  boundary  layer  also  contributes  losses  to  the  sound 
propagation  [Ref.  3,4].  While  these  viscous  and  thermal  effects  are  complicated,  if  the 
boundary  layer  thicknesses  are  small  compared  to  the  cross  sectional  dimensions  of  the 
resonator,  the  effects  can  be  included  by  a  modification  ofthe  wavenumber  k. 

The  frequency  dependent  viscous  and  thermal  boundary  layer  thickness  6^  and 


Sfr  are  given  by 


(2.83) 

(2.84) 
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where  //  is  the  viscosity,  K  is  the  thermal  conductivity,  c^is  the  specific  heat  and  p 
is  the  density  of  the  gas.  The  wave  number  k  is  then  given  by  [Ref.  5] 


(2.85) 


where  is  the  hydraulic  radius,  which  is  the  ratio  of  the  cross  sectional  area  to  the 
perimeter.  It  should  be  noted  that  since  the  hydraulic  radius  is  different  in  the  constriction 
and  in  the  main  channel,  and  that  k  is  now  complex.  The  imaginary  part  of  k  is  related 
to  the  thermo-viscous  damping. 


F.  END  EFFECTS 


In  the  earlier  analysis,  the  boundary  conditions  used  at  the  ends  of  the  constriction 
were  continuity  of  pressure  and  volume  velocity.  Using  only  plane  wave  analysis,  the 
pressure  and  particle  velocity  are  constant  across  the  cross  section.  However,  at  the 
constriction  end,  the  velocity  and  pressure  have  a  spatial  dependence  and  thus  the  plane 
wave  theory  is  not  as  accurate  near  the  ends  as  it  is  away  fi-om  them. 

The  spatial  variation  in  the  pressure  can  be  described  using  non-propagating 
higher  order  modes.  These  non-propagating  modes  store  energy  near  the  constriction 
ends.  Their  effect  on  the  plane  wave  propagation  can  be  modeled  as  a  lumped  impedance 
element.  [Ref  6],  The  lumped  impedance  can  be  obtained  fi'om  a  variety  of  methods. 
A  first  approximation  involves  modeling  the  impedance  as  an  added  mass  much  like 
the  end  correction  for  radiation  at  the  end  of  a  tube.  This  added  mass  can  be  computed 
using  conformal  mapping.  [Ref.  6].  A  more  accurate  method  involves  including  higher 
order  modes  near  the  constriction  end.  [Ref  7].  While  computationally  more  intensive 
than  the  conformal  mapping  solution,  the  higher  order  mode  solution  is  more  accurate. 

After  computing  the  lumped  impedance  Z(<o),  the  boundary  condition  for  the 
pressure  becomes  (the  volume  velocity  boundary  condition  is  unchanged.) 

P,(/)  =  P,(/)  +  Z(0)U2(/)  (2.86) 
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where  Pj  is  the  pressure  in  unconstricted  region  and  Pj  is  the  pressure  in  constriction. 
U2  is  the  velocity  in  the  constriction.  The  constriction  is  at  x  =  /. 
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m.  EXPERIMENTAL  APPARATUS  AND  PROCEDURE 


A.  INTRODUCTION 

A  description  of  the  experimental  apparatus  and  the  procedures  for  data 
acquisition  are  presented  in  this  chapter.  The  first  section  will  include  a  description  of  the 
annular  resonator,  constrictions,  the  driver,  the  microphone  and  electronic 
instrumentation.  The  second  section  discusses  the  data  acquisition  system.  The  third 
section  is  a  detailed  discussion  of  the  experimental  procedure. 

B.  EXPERIMENTAL  APPARATUS 

1.  Resonator  and  Driver 

Figure  3-1  is  an  illustration  of  the  annular  resonator  and  acoustic  driver.  The  sides 
and  bottom  of  resonator  are  made  of  aluminum.  The  top  of  resonator  is  made  of 
plexiglas  to  allow  visual  sighting  of  the  relative  positioning  of  the  constriction  and  the 
driver.  More  details  about  each  part  will  follow. 

The  dimensions  of  the  resonator  were  chosen  to  give  longitudinal  resonance 
fi-equencies  in  the  300 -500  Hz  range  and  to  be  well  below  the  cutoff  fi-equency 
for  cross  modes.  The  resonator  consists  of  two  main  parts  as  shown  on  Fig.  3-1. 
One  piece  is  an  aluminum  ring  with  15.30±0.02cm  inner  radius  and  1cm  thick  wall. 
This  ring  formed  the  outer  side  wall  of  the  resonator.  The  other  piece  is  a  solid 
aluminum  cylinder  which  is  concentric  with  the  outer  ring.  The  inner  radius  of 
resonator  is  10.0+0.01  cm.  The  height  of  resonator  is  5.0±0.01  cm.  The  resonator 
is  capped  on  the  open  side  with  1cm  thick  plexiglas  plate  with  a  radius  of  16.0  cm. 
The  bottom  plate  of  the  resonator  is  made  of  aluminum  and  its  thickness  is  1.0  cm. 
Sixteen  microphones  are  installed  in  the  cap  to  measure  the  spatial  acoustic  pressure 
distribution  inside  resonator.  The  installation  of  microphones  will  be  discussed  later. 
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The  cap  is  bolted  with  10.7  cm  long  threaded  rod  in  the  center  and  8  clamp  bars  as 
shown  on  Fig.  3-2.  To  ensure  a  tight  seal  (for  a  high  quality  factor),  vacuum  grease 
was  applied  to  the  surface  of  the  aluminum  ring  and  center  piece. 


The  acoustic  signal  was  generated  by  a  40  W  ,  16  ohm  external  driver.  The 
driver  was  connected  to  the  resonator  by  a  3.5  cm  long  piece  of  plastic  hose  with  0.4 
cm  diameter.  The  driver  was  excited  by  the  output  of  a  Hewlett  Packard  3 562A  dynamic 
signal  anal)^er  sent  through  a  Techron  5507  power  amplifier.  The  signal  analyzer 
was  used  to  determine  the  fi'equency  response  of  the  resonator.  Typically,  the 
resonator  was  driven  with  random  noise  with  a  1  Vrms  source  level.  To  measure  the 
mode  shape,  a  Hewlett  Packard  33 120 A  function  generator  was  used  to  excite  the 
resonator  at  the  resonance  fi'equency  determined  fi'om  the  built-in  pole-zero  analysis 
of  the  fi'equency  response  using  the  signal  analyzer. 
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Clamp  bar 


Threaded  Rod  &  nut 


Figure  3-2.  Illustration  of  clamp  bar,  threaded  rod  and  nut. 


2.  Microphones 

For  the  experiment,  16  Panasonic  WM60-AY  microphones  were  used  to 
measure  the  pressure  distribution  inside  resonator.  The  dimensions  of  the  microphones 
are  6  mm  in  height  and  5mm  in  diameter.  The  sensitivity  is  -42  dB  ±  3dB  re  IPa.  A 
Hewlett  Packard  623 7B  triple  output  power  supplier  provided  9  V  DC  to  each 
microphone. 


a.  Installation  of  microphones 

The  microphones  were  installed  as  shown  on  Fig.  3-3.  The  first 
intention  was  to  use  only  one  microphone  and  rotate  the  cap  for  each  different  position 
to  measure  the  mode  shape.  However,  the  rotation  of  cap  broke  the  sealing  between  the 
plexiglas  cap  and  aluminum  bodies  of  the  resonator  making  the  measurements  very  time 
consuming  and  nonrepeatable. 
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mic  1 


Figure  3-3.  Installation  of  16  microphones  on  the  plexiglas  cap. 

We  decided  to  install  16  microphones  in  order  to  protect  sealing  and  reduce  the  time 
required  for  data  acquisition.  For  the  precise  measurement  in  constricted  area,  6 
microphones  are  located  in  first  30  degree  region  where  a  constriction  is  located.  The 
remaining  10  microphones  are  located  30  degrees  apart.  The  microphones  were  installed 
in  the  center  of  the  annular  channel.  To  install  the  microphones,  16  holes  were  drilled 
part  way  through  the  plexiglas  cap.  These  holes  were  slightly  larger  than  the  diameter  of 
the  microphone.  A  smaller  hole  was  drilled  concentric  with  each  of  the  larger  holes.  This 
smaller  hole  penetrated  the  cap.  Therefore  the  hole  size  presented  to  the  inside  of  the 
resonator  was  much  smaller  than  the  microphone  diameter.  The  small  hole  diameter  is 
0.15  cm.  The  microphones  were  sealed  into  the  holes  with  silicon  glue.  The  error  in 
the  microphone  position  is  estimated  to  be  less  than  +  0.25°. 
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Built  in  speaker 
Chamber 


Figure  3-4 .  Chamber  used  for  microphone  calibration 

b.  Calibration  of  microphones 

Because  16  microphones  are  used  for  sensing  the  pressure  distribution,  the 
first  step  was  to  calibrate  all  the  microphones  over  the  300  ~  500  Hz  fi'equency  range 
interest.  For  the  calibration,  a  custom  made  chamber  with  a  built  in  speaker  was 
used.  As  shown  on  Fig.  3-4,  the  diameter  of  the  chamber  is  3.0  cm  and  the  height  is  2.50 
cm.  At  first  we  tried  to  measure  the  pressure  for  each  microphone  at  the  same 
position  relative  to  the  driver  in  an  otherwise  empty  resonator.  However,  the  phase 
shift  caused  by  the  resonators’  resonance  made  calibration  very  sensitive  to  small 
errors  in  position.  Instead,  the  top  cap  has  removed  and  the  calibration  chamber 
was  positioned  directly  over  each  microphone.  The  frequency  response  was  measured 
for  each  microphone  over  a  300  ~  500  Hz  range.  A  t5q)ical  fi'equency  response  is 
shown  in  Fig.  3-5.  This  sweeping  and  data  storage  were  done  by  Hewlett  Packard 
35665A  dynamic  signal  analyzer.  For  the  calibration,  the  signal  analyzer  was  set  to  swept 
sine  mode  with  a  settling  time  of  20  cycles  and  integration  time  of  50  cycles.  The 
chamber  was  placed  onto  the  cap  where  each  microphone  is  located  and  the  chamber 
was  sealed  on  the  cap  with  vacuum  grease. 
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O  X;353  Hz  Y;474.0161  m 


Figure  3-5.  Frequency  response  of  Nfic.  #1  covered  by  the  cahbration  chamber. 

Microphone#!  was  used  as  the  reference  for  comparison  calibrations  of  the 
other  15  microphones.  Each  microphone’s  output  was  divided  by  that  of 
microphone  #1  to  create  a  correction  factor  as  shown  on  Fig.  3-6.  When  measuring 
the  mode  shape,  the  output  from  each  microphone  was  divided  by  the  appropriate 
correction  factor  to  obtain  the  corrected  amplitude  at  a  particular  location  and 
driving  frequency.  For  example,  consider  the  case  when  the  measured  amplitude  from 
microphone  #  12  is  50  mV  at  455  Hz.  The  correction  factor  is  0.739  from  the  Fig. 
3-6.  Dividing  50  mV  by  0.739  gives  67.66  mV  as  the  corrected  measurement. 
Generally,  the  correction  factors  were  between  0.7  and  1.2. 

OX:455  Hz  Y:738.7183m 


300  Hz  500  Hz 

Figure  3-6.  Plot  of  correction  factor  for  Mic.  #12  referenced  to  Mic.  #l.The  correction 


factor  is  0.739  at  455  Hz. 


3. 


Constriction 


To  explore  a  wide  range  of  parameters,  various  sizes  of  constriction  are  necessaiy. 
When  considering  the  material  selection  for  this  constriction,  a  couple  of  factors  came 
to  mind.  The  material  had  to  be  rigid  and  have  good  machining  propertied  for  ease 
of  fabrication.  These  two  factors  led  to  the  selection  of  PVC  (Polyvinylchloride). 
Figure  3-7  shows  the  design  along  with  the  nominal  and  measured  dimensions  for  the 
various  constrictions  used.  For  the  experiment,  3  different  angular  widths  were  used. 
The  width  of  constriction  is  given  by  the  0  shown  in  Fig.  3-7.  The  three  widths 

are  22.5°  ,  11.25°  and  5.625°.  The  22.5°  constriction  corresponds  to  of  the 

16 

effective  circumference  of  the  resonator.  11.25°  corresponds  to  of  and 

5.625°  corresponds  to  —  of  The  area  ratio  is  the  ratio  of  open  area  in  the 

64 

constriction  (which  is  the  shaded  black  square  in  Fig.  3-7)  to  cross  sectional  area  of 
resonator.  In  Fig  3-7,  SI  stands  for  the  open  area  and  S2  stands  for  the  open  cross 
sectional  area  of  the  constriction.  The  table  in  Fig.  3-7  gives  the  geometrical  parameters 
of  32  different  constrictions. 
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Figure  3-7a.  Nominal  geometry  of  the  constrictions 


unit:  cm 


e 

cn^ ) 

area  ratio(Sl/S2) 

22.53  ±  0.06 

2.5911 

26.447 

0.0980 

5.2441 

0.1983 

7.8400 

0.2964 

13.1769 

0.4982 

18.0170 

0.6812 

11.30±  0.05 

2.5911 

26.447 

0.0980 

5.2441 

0.1983 

7.8400 

0.2964 

13.1769 

0.4982 

18.0170 

0.6812 

5.65  ±  0.02 

2.5911 

26.447 

0.0980 

5.2441 

0.1983 

7.8400 

0.2964 

13.1769 

0.4982 

18.0170 

0.6812 

Figure  3-7b.  Measured  geometry  of  the  constrictions. 


4.  Electronic  Instrumentation 


Figure  3-8  is  a  block  diagram  of  the  electronic  instrumentation  used  during 
the  measurements.  A  Hewlett  Packard  3 562 A  Dynamic  signal  analyzer  is  used  to 
determine  the  frequency  response  in  the  band  from  300  Hz  to  500  Hz.  The  signal 
from  this  analyzer  is  an^lified  by  a  Techron  5507  power  supply  amplifier  and  used  as 
the  input  signal  to  the  driver.  The  pole-zero  analysis  function  built  into  the  signal 
analyzer  is  used  to  extract  the  resonance  frequencies  from  the  measured  frequency 
response.  Once  the  resonance  frequency  is  found,  a  Hewlett  Packard  33120A 
Function/Arbitrary  Waveform  Generator  is  used  to  generate  a  1  Vpp  sinusoidal 
signal  at  the  obtained  resonance  frequency.  The  signal  is  also  an:5)lified  by  the 
same  Techron  anplifier  and  used  as  the  input  signal  to  the  driver.  Either  the  signal 
from  the  analyzer  or  the  function  generator  signal  is  monitored  on  a  Tektronics  2336 
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Oscilloscope.  The  outputs  from  the  microphones  are  monitored  on  the  same  oscilloscope. 
A  custom  made  microphone  switch  board  selects  each  microphone  to  measure  the 
pressure  anplitude  at  different  positions.  The  microphone  output  anq)litude  is  measured 
using  a  Standford  Research  530  Lock  Li  an^lifier.  In  order  to  measure  the 
tenqierature  of  the  resonator,  an  Omega  450  AET  thermocouple  thermometer  type  E  is 
used.  The  thermocouple  was  calibrated  to  ±0.1  ®C  using  a  distilled  water  ice  bath.  The 
thermocouple  is  positioned  on  the  surface  of  the  outer  aluminum  ring  of  resonator. 
The  tenqierature  inside  resonator  is  measured  frequently  and  be  recorded  whenever  there 
is  variation. 
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Lock  in  amplifier 
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Figure  3-8.  Block  diagram  of  instrumentation  setup. 


C.  EXPERIMENTAL  PROCEDURE 


1.  Resonator  Setup 

Prior  to  any  data  acquisition,  the  resonator  is  assembled  with  the  constriction  in 
the  desired  position  relative  to  the  driver.  For  the  experiment,  three  major  positions  are 
used.  They  are  45®,  90°  and  180°  from  driver.  The  position  is  measured  from  the 
driver  to  the  center  of  constriction.  The  top  cap  is  oriented  so  that  0°  is  set  at  the 
left  edge  of  constriction.  The  0°  setting  is  performed  by  rotation  of  the  cap  by  hand. 
For  the  error  measurement  of  0°  setting,  three  people  conducted  rotation  of  the 
cap.  Each  person  tried  three  times  and  we  found  that  the  error  by  rotation  was  less 
than  +0.5°.  So,  with  the  mechanical  error  ±0.25°  in  drilling  the  microphone  hole,  the 
total  error  of  microphone  positioning  is  less  than  +0.75°.  With  the  constriction 
positioned,  the  proper  amount  of  vacuum  grease  is  put  on  the  open  surface  of  resonator 
and  the  top  cap  is  secured  with  the  clamp  bars. 

2.  Determining  Resonator  Characteristics 

The  Hewlett  Packard  3 562 A  Dynamic  Signal  Analyzer  is  used  to  characterize 
the  resonator.  The  swept  sine  mode  of  the  analyzer  provides  a  sine  wave  signal  to 
the  system  at  a  given  range  of  frequencies.  The  frequency  response  of  the  system  is 
used  to  determine  the  resonance  frequency  and  quality  factor  of  several  modes  of  the 
resonator.  When  only  one  mode  is  excited,  the  peak  pressure  amplitude  frequency  is 
same  as  the  resonance  frequency  given  by  pole-zero  analysis.  When  two  modes  are 
excited,  pole-zero  analysis  is  used  to  determine  the  resonance  frequency.  The  signal 
analyzer  gives  poles  in  the  fbrm  a  +  jb  .  The  resonance  frequency  is  simply  the 

value  of  b .  The  quality  factor  can  be  also  obtained  by  ^  from  that  form.  This 

method  has  been  explained  in  the  theory  chapter.  These  measurements  are  performed 
for  each  constriction  type  and  each  relative  position  of  constriction  and  driver. 
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3. 


Data  Collection 


Data  collection  is  performed  in  sets  of  the  positions  and  different  sizes  of 
constriction.  Before  data  collection,  the  initial  ambient  temperature  and  noise  level 
is  checked.  The  output  of  the  Hewlett  Packard  6237B  Triple  output  power 
supply  is  also  set  to  DC  9V.  The  measurement  is  then  commenced  by  getting  the 
resonance  frequency  from  the  Hewlett  Packard  3 562 A  Dynamic  signal  analyzer.  The 
microphone  switch  board  is  set  to  the  microphone  which  is  closest  to  the  driver.  Then 
the  resonance  frequency  and  the  temperature  of  resonator  was  measured.  The 
procedure  to  acquire  resonance  frequency  has  been  mentioned  in  subsection  2  in 
this  section.  Getting  quality  factor  also  follows  the  same  procedure  previously 
mentioned.  After  finding  resonance  frequency,  a  data  run  is  commenced  by 
hooking  up  the  Hewlett  Packard  33120A  fimction  generator  to  the  input  amplifier  for 
driver  and  setting  the  frequency  and  voltage.  The  driving  frequency  should  be  the 
resonance  frequency  obtained  by  previous  procedure  and  the  voltage  is  1  Vpp.  The 
acoustic  amplitude  is  measured  with  each  microphone.  With  the  driver  still  running,  the 
switch  board  selects  a  particular  microphone  and  the  value  of  amplitude  at  each 
position  is  measured  with  the  SR  530  Lock  -  In  amplifier.  At  this  time,  if  there  is 
variation  in  the  resonator  temperature,  it  is  noted.  All  data  recorded  from  the  16 
microphones  are  scaled  by  the  correction  factor  obtained  from  microphone  calibration. 
The  corrected  data  go  to  computer  program  to  plot  mode  shape.  This  concludes  data 
acquisition  of  one  mode. 
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Procedure  for  Wave  Guide  Characterization  Collection 


1.  Date:  /  /97  (  ~  ) 

2.  Ambient  Condition 

Room  Temper :  °C  (  K)  /  Speed  of  Sound  =  m/s 

Ambient  Pressure:  mbar  /  Initial  Temp  of  resonator:  ®C 


3.  Configuration  (Constriction) 

#  of  Const.  :  EA  /  Size : 

Area  Ratio : 

Position  of  Const. :  ®  (CCW/CW)  from  driver 

Position  of  Driver  :  ® 


4.  Resonance  Freq: 


/  Q:  /Mode: 


5.  Source  information 

Input :  Vpp 


/  Type  of  source  : 


6.  Data  :WG  with _  Const(s)  =>  T _  (  mode) 


Receiver 

Vrms  of 
receiver 

Corerected 

Data 

Temp  of 
Torus 

Resonance 

Freq 

Q 

1  (000) 

( 

) 

2  (006) 

( 

) 

3  (012) 

( 

) 

4  (018) 

( 

) 

5  (024) 

( 

) 

6  (030) 

( 

) 

7  (060) 

( 

) 

8  (090) 

( 

) 

9  (120) 

( 

) 

10  (150) 

( 

) 

11  (180) 

( 

) 

12  (210) 

( 

) 

13  (240) 

( 

) 

14  (270) 

( 

) 

15  (300) 

( 

) 

16  (330) 

( 

) 

Figure  3-9.  Data  collection  form. 
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IV.  RESULTS  AND  DISCUSSION 


The  purpose  of  this  chapter  is  to  show  the  comparison  between  theory  and 
experimental  results.  The  result  of  resonance  frequency  measurements  are  discussed  first 
and  then  the  mode  shapes. 


A.  RESONANCE  FREQUENCIES 


The  measured  resonance  frequencies  are  determined  from  pole-zero  analysis  of  the 
frequency  response.  Because  the  resonance  frequencies  were  measured  at  different 
temperatures,  they  were  corrected  to  20®C  in  order  to  make  comparison  of  data 
recorded  at  various  temperatures  with  the  theoretical  resonance  frequences  more 
straightforward.  The  equation  for  conversion  is 


f -20  ~  ft 


measured 


1  + 


20 

273 


1  + 


^measured 

273 


(4.1) 


The  theoretical  resonance  frequencies  were  calculated  for  two  different  cases.  One  case 
was  without  end  corrections  and  the  other  was  with  end  corrections.  All  the  results  are 
listed  in  Tables  4-1  and  4-2.  Tables  4-la  and  4.2a  list  the  resonance  frequencies  for  the 
low  (L)  and  high  (H)  modes  for  the  various  constriction  lengths  and  area  ratios.  Tables  4- 
Ib  and  4-2b  lists  the  temperatures  at  which  the  frequencies  were  measured.  Tables  4-lc 
and  4-2c  list  the  20°C-corrected  frequencies.  Tables  a,  b  and  c  are  identical  for  Tables  4-1 
and  4-2.  Tables  4-ld  and  4-2d  give  the  calculated  resonance  frequencies.  Table  4-ld 
ignores  end  corrections.  Table  4-2d  includes  end  corrections  as  calculated  from  higher 
order  non  propagating  modes.  Tables  4-le  and  4-2e  list  the  percent  error  between 
measured  and  calculated  resonance  frequencies. 

It  is  evident  that  the  inclusion  of  end  corrections  has  little  effect  on  the  prediction 
for  the  high  modes.  In  both  cases  the  average  error  is  approximately  0.15%.  However, 
end  corrections  are  very  important  for  the  low  mode.  This  point  is  evident  from  Fig.  4-1 
where  it  is  seen  that  inclusion  of  end  corrections  is  important  for  the  low  mode  only. 
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Figure.  4-1  shows  the  20®C-corrected  resonance  frequencies  as  well  the  predicted  results 
(without  and  with  end  corrections)  as  functions  of  area  ratio. 

In  this  section  it  was  shown  that  end  corrections  are  an  important  factor  to  include 
when  calumniating  resonance  frequencies,  especially  for  the  low  mode.  In  the  next  section, 
it  is  shown  that  they  are  also  important  when  calculating  mode  shapes. 
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1.  Comparison  of  Measured  Resonance  Frequency  with  Calculated  Results 
without  End  Corrections 


Table  4- la.  Measured  resonance  frequencies  for  different  constrictions. 


AR  0.1  1  0.2  0.3  0.5  0.7  1.0 


H  L  I  H  L  H  L  I  H  L  H 


L/16 

22.5 

22.5 

21.4 

21.4 

21.8 

21.8 

22.0 

22.0 

22.4 

22.4 

22.1 

L/32 

23.9 

23.9 

23.3 

23.3 

23.6 

23.6 

23.7 

23.7 

21.6 

21.6 

L/64 

24.0 

24.0 

22.3 

22.3 

22.3 

22.3 

22.3 

22.3 

23.5 

23.5 

Table  4-lb.  Temperatures  for  each  measurement  from  table  4-la. 


Table  4-lc.  Measured  resonance  frequencies  corrected  to  20°C. 


Table  4- Id.  Calculated  resonance  frequencies  at  20®C,  without  end  corrections. 


AR  0.1  0.2  0.3  0.5  0.7  1.0 

si»  L  I  H  L  I  H  L  r~a  L  I  H  ~L  I  H  _ 


L/16 

3.48 

0.07 

3.71 

0.18 

2.91 

0.21 

1.29 

0.13 

0.32 

0.04 

0.05 

L/32 

6.98 

0.16 

5.65 

0.19 

3.83 

0.14 

1.53 

0.13 

0.33 

0.10 

L/64 

10.1 

0.13 

6.54 

0.20 

4.19 

0.20 

1.44 

0.20 

0.38 

0.11 

Table  4-le.  Percent  error  between  measured  and  calculated  resonance  frequency. 


45 


2.  Comparison  of  Measured  Resonance  Frequency  with  Calculated  Results  with 
End  Corrections 


Table  4-2a.  Measured  resonance  frequencies  for  different  constrictions. 


Table  4-2b.  Temperatures  for  each  measurement  from  table  4-2a. 


Table  4-2c.  Measured  resonance  frequencies  corrected  to  20°C. 


Table  4-2d.  Calculated  resonance  frequencies  at  20®C  with  end  corrections. 


Table  4-2e.  Percent  error  between  measured  and  calculated  resonance  frequency. 
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measured  &  calculated  resoriance  freq  vs  area  ratio 


Figure  4- la.  Comparison  of  calculated  and  measured  resonance  frequencies  versus  area 
ratio  for  the  L/16  constriction. 


measured  &  calculated  resonance  freq  vs  area  ratio 


Figure  4-lb.  Comparison  of  calculated  and  measured  resonance  frequencies  versus  area 
ratio  for  the  L/32  constriction. 
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measured  &  calculated  resonance  freq  vs  area  ratio 


Figure  4-lc.  Comparison  of  calculated  and  measured  resonance  frequencies  versus  area 
ratio  for  the  L/64  constriction. 

B.  MODE  SHAPES 

Figures  4-2a  and  b  show  the  comparison  of  the  measured  and  computed  mode 
shapes  for  the  L/16,  0.1  area  ratio  constriction.  The  measurements  have  been  normalized 
to  the  highest  value  of  pressure  amplitude.  The  amplitude  of  the  mode  shape  is 
determined  from  least  squares  fit  to  the  measured  data.  The  driver  is  located  45®  from  the 
constriction.  Figure  4-2a  shows  the  results  for  the  high  mode.  It  is  seen  that  there  is 
good  agreement  between  the  measured  and  predicted  mode  shapes,  even  though  end 
corrections  are  not  included  in  the  calculation.  The  driver  is  located  at  a  place  where  both 
modes  should  be  excited.  However,  the  modes  are  well  separated  with  this  constriction, 
so  there  is  little  excitation  of  the  low  mode.  Single  mode  excitation  is  indicated  by  a  high 
standing  wave  ratio  from  the  mode  shape. 

Figure  4-2b  shows  the  mode  shape  for  the  low  mode.  The  agreement  between  the 
measured  and  calculated  mode  shape  is  significantly  worse  than  with  the  high  mode. 
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Figures  4-2c  and  d  show  the  mode  shapes  agmn  but  this  time  end  corrections  are 
taken  into  account.  Comparing  Fig.  4-2c  and  4-2a,  it  is  seen  that  end  corrections  have 
little  effect  on  the  high  mode.  However,  they  are  very  important  for  the  low  mode,  just 
as  was  discussed  in  the  previous  section.  Figure  4-2d  shows  that  when  end  corrections 
are  taken  into  account,  the  agreement  between  calculated  and  measured  mode  shape 
improves  significantly. 

Figures  4-3a  and  b  show  the  mode  shapes  for  the  case  when  the  driver  is  located 
45®  from  the  L/32,  0.7  area  ratio  constriction.  In  the  case  the  constriction  is  relatively 
short  and  the  area  ratio  is  relatively  large.  As  a  result,  the  modes  are  not  well  separated. 
So  although  the  driver  is  exdted  at  the  resonance  frequency  of  one  mode,  there  is 
significant  excitation  of  the  other  mode,  as  evidenced  by  the  low  standing  wave  ratio. 

Figure  4-3a  shows  the  mode  shape  when  the  resonator  is  driven  at  the  high  mode 
resonance  frequency.  End  correction  are  ignored  in  the  calculation.  The  agreement 
between  measured  and  calculated  mode  shape  is  not  as  good  as  in  the  previous,  well 
separated-mode  case  (Fig.  4-2).  Figure  4-3b  shows  the  results  for  the  low  mode.  Here 
the  agreement  between  theory  and  experiment  is  worse.  In  both  cases,  the  standing  wave 
ratio  is  noticeably  lower  that  in  the  well-separated  case. 

Figure  4-3  c  and  d  show  the  effect  of  including  end  corrections.  There  is  some 
improvement  in  the  agreement  for  the  high  mode  case.  There  is  considerable  improvement 
in  the  agreement  for  the  low  mode,  and  some  improvement  for  the  high  mode.  Many  other 
comparisons  are  provided  in  the  Appendix.  The  cases  presented  in  this  chapter  are 
representative.  The  overall  agreement  in  mode  shape  is  good  even  in  the  overlapping¬ 
mode  case. 

The  results  of  this  chapter  indicate  that  end  corrections  are  important  both  for 
calculation  of  resonance  frequency  and  mode  shapes.  Also  end  corrections  are  more 
important  for  the  low  mode.  This  is  because  the  acoustic  impedance  of  the  low  mode  is 
low  in  comparison  to  that  of  the  high  mode.  The  addition  of  the  effective  impedance  of 
the  end  corrections  is  of  more  consequence,  then,  to  the  low  mode. 
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Mode  Shape  for  Length  L/16,  Area  Ratio  0.1,  and  f=461.5 


Figure  4-2a.  Mode  shape  for  the  high  mode  when  the  driver  is  located  45°  from  the 
constriction  having  length  of  L/16  and  area  ratio  of  0.1.  End  corrections  are  not  included 
in  the  calculation. 


Mode  Shape  for  Length  L/16,  Area  Ratio  0.1,  and  f=308.2 


Figure  4-2b.  Mode  shape  for  the  low  mode  when  the  driver  is  located  45°  from  the 
constriction  having  length  of  L/16  and  area  ratio  of  0.1.  End  corrections  are  not  included 


in  the  calculation. 


Mode  Shape  for  Length  L/16,  Area  Ratio  0.1,  and  f=461.5 


Figure  4-2c.  Mode  shape  for  the  high  mode  when  the  driver  is  located  45°  from  the 
constriction  ha\dng  length  of  L/16  and  area  ratio  of  0.1.  End  corrections  are  included  in 
the  calculation. 


Mode  Shape  for  Length  L/16.  Area  Ratio  0.1,  and  f=308.2 


Figure  4-2d.  Mode  shape  for  the  low  mode  when  the  driver  is  located  45°  from  the 
constriction  having  length  of  L/16  and  area  ratio  of  0.1.  End  corrections  are  included. 
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Mode  Shape  for  Length  L/32,  Area  Ratio  0.7,  and  f=439.3 


Figure  4-3a.  Mode  shape  for  the  high  mode  when  the  driver  is  located  45®  from  the 
constriction  having  length  of  L/32  and  area  ratio  of  0.7.  End  corrections  are  not  included. 


Mode  Shape  for  Length  L/32,  Area  Ratio  0.7,  and  f=426.9 


Figure  4-3b.  Mode  shape  for  the  low  mode  when  the  driver  is  located  45°  from  the 
constriction  having  length  of  L/32  and  area  ratio  of  0.7.  End  corrections  are  not  included. 
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Mode  Shape  for  Length  L/32,  Area  Ratio  0.7,  and  f=439.3 


Figure  4-3c.  Mode  shape  for  the  high  mode  when  the  driver  is  located  45°  from  the 
constriction  having  length  of  L/32  and  area  ratio  of  0.7.  End  corrections  are  included 


Mode  Shape  for  Length  L/32,  Area  Ratio  0.7,  and  f=426.9 


Figure  4-3d.  Mode  shape  for  the  low  mode  when  the  driver  is  located  45°  from  the 
constriction  having  length  of  L/32  and  area  ratio  of  0.7.  End  corrections  are  included 
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V. 


SUMMARY  AND  CONCLUSIONS 


The  purpose  of  this  thesis  is  to  investigate  frequency  splitting  and  mode  shapes  in  a 
constricted  annular  resonator.  The  measurements  are  intended  to  serve  as  preliminary 
verification  of  a  model  predicting  the  performance  of  an  annular  thermoacoustic  prime 
mover. 

The  constriction  splits  the  degeneracy  of  the  plane  wave-like  modes  of  a  uniform 
cross  section  annular  resonator.  One  mode  has  a  pressure  antinode  at  the  constriction,  the 
other  a  pressure  node.  The  resonance  frequency  of  the  former  is  higher  than  that  of  the 
latter.  The  difference  in  resonance  frequency  depends  upon  the  length  and  porosity  of  the 
constriction.  The  effects  of  dissipation  cause  the  frequency  responses  of  the  individual 
modes  to  overlap. 

An  equivalent  plane  wave  model  is  developed  to  predict  the  resonance  frequencies 
and  mode  shapes.  The  experimental  resonance  frequencies  are  determined  from  pole-zero 
analysis  of  the  measured  frequency  response. 

It  is  shown  that  it  is  critical  to  take  into  account  end  corrections  at  the  constriction 
to  get  good  agreement  between  the  measured  and  predicted  results.  End  corrections  are 
more  important  for  the  low  mode  than  the  high  mode. 

In  general  the  percent  error  between  measured  and  predicted  resonance 
frequencies  is  less  than  0.2  %  for  high  mode  and  0.3  %  for  the  low  mode.  The  measured 
mode  shape  are  also  in  good  agreement  with  predictions. 

The  results  of  this  thesis  indicate  that  the  simple  constricted  annular  resonator  is 
well  understood.  The  sensitivity  of  the  resonance  frequencies  and  mode  shapes  to  end 
corrections  indicates  that  the  annular  resonator  is  a  good  geometry  for  investigating  these 
effects  further.  For  instance,  how  does  one  model  the  end  corrections  for  a  thermoacoustic 
stack?  Are  the  end  corrections  for  a  stack  and  a  constriction  having  the  same  porosity 
substantially  different?  These  questions  are  left  for  future  work. 
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APPENDIX  DATA  TABLE  &  MODE  SHAPE 


Constriction 


length;  L/ 16  ,  area  ratio;  0.1  ,  angle  from  driver;  180° 


463.7  Hz 


Q;109 


1  (000) 


2  (006) 


3  (012) 


4  (018) 


5  (024) 


6  (030) 


7  (060) 


8  (090) 


Vrms 

Corrected 

(mV) 

Data 

70.5 

70.5 

67.3 

69.7 

55.9 

71.5 

52.8 

70.5 

44.3 

69.4 

58.9 

69.1 

39.2 

53.1 

15.9 

20.8 

Mode;  high 


Vrms 

(mV) 


13.1 


33.6 


79.4 


45.4 


33.1 


5.8 


27.4 


48.6 


9  (120) 


10  (150) 


11  (180) 


12  (210) 


13  (240) 


14  (270) 


15  (300) 


16  (330) 


Corrected 

Data 


16.6 


50.1 


67.0 


65.1 


43.0 


Mode  Shape  for  Length  L/16,  Area  Ratio  0.1,  and  f=463.7 


150  200  25 

Theta  (degrees) 


Experiment  2. 


Experiment  3. 


Constriction 


length:  L/ 16  ,  area  ratio:  0.1  ,  angle  from  driver:4y 


461.5  Hz 


1  (000) 


2  (006) 


3  (012) 


4  (018) 


5  (024) 


6  (030) 


7  (060) 


8  (090) 


Vrms 

Corrected 

(mV) 

Data 

49.8 

49.8 

48.7 

50.4 

39.8 

50.9 

37.8 

50.5 

31.4 

49.2 

41.7 

48.9 

27.8 

35.5 

11.3 

14.8 

Mode:  high 


Vrms 

(mV) 


9  (120) 


10  (150) 


11  (180) 


12  (210) 


13  (240) 


14  (270) 


15  (300) 


16  (330) 


Corrected 

Data 


11.9 


36.1 


48.5 


46.8 


30.8 


5.5 


20.7 


41.4 


Mode  Shape  for  Length  L/16,  Area  Ratio  0.1,  and  f=461.5 


— Theory 
^  Measurement 
Edge  of  Const. 
Driver  Position 


Theta  (degrees) 


Experiment  4. 


Experiment  5. 


Experiment  6. 


Constriction 


Resonance  frq. 


1  (000) 


2  (006) 


3  (012) 


4  (018) 


5  (024) 


6  (030) 


7  (060) 


8  (090) 


length:!/ 32  ,  area  ratio:  0.1  ,  angle  from  driver:90® 


328.7  Hz 


Vrms 

(mV) 


20.3 


3.8 


21.2 


27.9 


26.1 


37.4 


36.7 


35.6 


Q:  69.4 


Corrected 

Data 


20.3 


3.9 


26.7 


36.4 


40.8 


43.6 


49.7 


45.9 


Mode:  low 


9  (120) 


10  (150) 


11  (180) 


12  (210) 


13  (240) 


14  (270) 


15  (300) 


16  (330) 


Vrms 

(mV) 

Corrected 

Data 

28.7 

35.9 

15.0 

22.2 

4.0 

3.4 

10.9 

15.5 

24.8 

32.2 

31.6 

42.3 

46.0 

48.3 

38.9 

45.8 

Mode  Shape  for  Length  L/32,  Area  Ratio  0.1 ,  and  f=328.7 


— Theory 
^  Measurement 
Edge  of  Const. 
Driver  Position 


Theta  (degrees) 


Experiment  7. 


Experiment  8. 


Mode  Shape  for  Length  L/32,  Area  Ratio  0.1 ,  and  f=328.7 


Experiment  9. 


Constriction 

length:  L/64  ,  area  ratio:  0.1  ,  angle  from  driver:  180° 

Resonance  frq.:  443.2  Hz 

Nficrophone 

1 m 

Vrms 

(mV) 

Corrected 

Data 

Microphone 
#  (0) 

Vrms 

(mV) 

Corrected 

Data 

1  (000) 

72.8 

72.8 

9  (120) 

24.9 

31.4 

2  (006) 

67.7 

69.9 

10  (150) 

41.4 

61.7 

3  (012) 

57.5 

73.3 

11  (180) 

85.9 

72.6 

4  (018) 

53.6 

71.5 

12  (210) 

45.1 

64.7 

5  (024) 

43.9 

68.7 

13  (240) 

29.7 

38.5 

6  (030) 

56.9 

66.6 

14  (270) 

1.4 

1.9 

7  (060) 

31.1 

42.0 

15  (300) 

33.0 

34.1 

8  (090) 

4.0 

5.2 

16  (330) 

52.9 

62.2 

Mode  Shape  for  Length  L/64,  Area  Ratio  0.1,  and  f=443.2 
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Experiment  10. 


Constriction  length:  1 

[764  ,  area  ratio:  0. 1  ,  angle  from  driver: 90® 

Resonance  frq.:  350.2  Hz 

Q;  41  Mode:  1 

Nficrophone 
#  O) 

Vrms 

(mV) 

Corrected 

Data 

Nficrophone 
#  (0) 

Vrms 

(mV) 

Corrected 

Data 

1  (000) 

7.1 

7.1 

9  (120) 

17.9 

22.4 

2  (006) 

8.5 

8.7 

10  (150) 

9.2 

13.6 

3  (012) 

13.8 

17.4 

11  (180) 

2.0 

1.7 

4  (018) 

15.8 

20.7 

12  (210) 

7.3 

10.4 

5  (024) 

14.6 

22.8 

13  (240) 

16.0 

20.8 

6  (030) 

21.2 

24.8 

14  (270) 

19.8 

26.4 

7  (060) 

21.9 

29.5 

15  (300) 

28.0 

29.3 

8  (090) 

21.8 

28.2 

16  (330) 

22.3 

26.2 

Mode  Shape  for  Length  L/64,  Area  Ratio  0.1 ,  and  f=350.2 
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Constriction 


Resonance  frq. 


length;  L/ 64  ,  area  ratio:  0. 1  ,  angle  from  driver:45' 


1  (000) 


2  (006) 


3  (012) 


4  (018) 


5  (024) 


6  (030) 


7  (060) 


8  (090) 


443.1Hz 


Vrms 

(mV) 


39.8 


37.0 


31.1 


29.1 


23.9 


31.0 


16.8 


2.2 


Q:  87 


Corrected 

Data 


39.8 


38.2 


39.7 


38.8 


37.4 


36.3. 


22.7 


2.9 


9  (120) 


10  (150) 


11  (180) 


12  (210) 


13  (240) 


14  (270) 


15  (300) 


16  (330) 


Vrms 

(mV) 

Corrected 

Data 

13.9 

17.5 

22.7 

33.8 

47.0 

39.7 

24.7 

35.4 

16.1 

20.9 

0.8 

1.1 

18.6 

19.2 

29.2 

34.3 

Mode  Shape  for  Length  L/64,  Area  Ratio  0.1,  and  f=443.1 


100  150 

Theta  (degrees) 


Experiment  12. 


Mode  Shape  for  Length  L/64,  Area  Ratio  0.1 ,  and  f-348.6 


Experiment  13. 


Constriction 


length:  L/ 16  ,  area  ratio:  0.7  ,  angle  from  driver:  1 80° 


444.2  Hz 


1  (000) 


2  (006) 


3  (012) 


4  (018) 


5  (024) 


6  (030) 


7  (060) 


8  (090) 


Vrms 

Corrected 

(mV) 

Data 

69.3 

69.3 

67.6 

69.9 

55.9 

71.3 

53.1 

70.8 

44.5 

69.6 

58.7 

68.7 

36.8 

49.8 

13.0 

17.0 

Mode:  high 


50  100  150  200  250  300  350 

Theta  (degrees) 
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Experiment  14. 


Mode  Shape  for  Length  L/16,  Area  Ratio  0.7,  and  f-421.2 


Experiment  15. 


Constriction 


Resonance  frq. 


length:  L/  16  ,  area  ratio:  0.7  ,  angle  from  driver:45‘ 


1 

(000) 

2 

(006) 

3 

(012) 

4 

(018) 

5 

(024) 

6 

(030) 

7 

(060) 

8 

(090) 

444.0  Hz 


Vrms 

(mV) 


41.6 


40.5 


33.1 


31.5 


26.3 


34.5 


21.4 


.5 


Q:93 


Corrected 

Data 


41.6 


41.9 


42.2 


42.0 


41.2 


40.4 


29.0 


9.8 


Mode:  high 


9  (120) 


10  (150) 


11  (180) 


12  (210) 


13  (240) 


14  (270) 


15  (300) 


16  (330) 


Vrms 

(mV) 

Corrected 

Data 

10.7 

13.5 

21.5 

32.0 

48.3 

40.8 

27.6 

39.6 

20.6 

26.7 

IBB 

BIB 

15.1 

15.6 

28.2 

33.1 

Mode  Shape  for  Length  L/16,  Area  Ratio  0.7,  and  f=444.0 


— Theory 
^  Measurement 
Edge  of  Const. 
Driver  Position 


0  50  100  150  200  250  3 

Theta  (degrees) 


I  ...  -m 


Experiment  16. 


I  Constriction  length:  L/ 16  ,  area  ratio:  0.7  ,  angle  from  (iriver:45° 

1  Resonance  frq.:  421.9  Hz 

Microphone 
#  O) 

Vrms 

(mV) 

Corrected 

Data 

Microphone 
#  (9) 

Vrms 

(mV) 

Corrected 

Data 

1  (000) 

14.3 

14.3 

9  (120) 

43.4 

54.5 

2  (006) 

6.7 

6.9 

10  (150) 

25.8 

38.4 

3  (012) 

3.1 

3.9 

11  (180) 

12.7 

10.8 

4  (018) 

8.4 

11.2 

12  (210) 

14.4 

20.6 

5  (024) 

11.9 

18.6 

13  (240) 

34.9 

45.3 

6  (030) 

21.9 

25.6 

14  (270) 

42.4 

56.2 

7  (060) 

36.7 

49.5 

15  (300) 

54.6 

56.5 

8  (090) 

45.9 

59.8 

16  (330) 

35.7 

41.9 

Mode  Shape  for  Length  L/16,  Area  Ratio  0.7,  and  f=421.9 
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Experiment  17. 


Constriction 


length:  L/  32  ,  area  ratio:  0.7  ,  angle  from  driver:  180® 


441.1Hz 


Q:  114 


1  (000) 


2  (006) 


3  (012) 


4  (018) 


5  (024) 


6  (030) 


7  (060) 


8  (090) 


Vrms 

(mV) 

Corrected 

Data 

69.4 

69.4 

67.4 

69.6 

55.5 

70.8 

51.5 

68.7 

43.6 

68.2 

56.8 

66.5 

32.8 

44.3 

7.8 

10.2 

Mode:  high 


Vrms 

(mV) 


20.4 


37.5 


81.4 


45.5 


32.4 


9  (120) 


10  (150) 


11  (180) 


12  (210) 


13  (240) 


14  (270) 


15  (300) 


16  (330) 


Corrected 

Data 


25.7 


55.9 


68.8 


65.2 


42.0 


.6 


27.6 


56.4 


Mode  Shape  for  Length  L/32,  Area  Ratio  0.7,  and  ^=441.1 


Experiment  18. 


Constriction  length:  ] 

U  32  ,  area  ratio:  0.7  ,  anj 

?le  from  diiver:90® 

Resonance  frq.:  426.2  Hz 

Q:  129  Mode:  1 

Microphone 
#  (9) 

Vrms 

(mV) 

Corrected 

Data 

Mcrophone 

#  (9) 

Vrms 

(mV) 

Corrected 

Data 

1  (000) 

8.7 

8.7 

9  (120) 

57.8 

72.6 

2  (006) 

2.8 

2.9 

10  (150) 

31.7 

47.2 

3  (012) 

10.7 

13.6 

11  (180) 

7.2 

6.1 

4  (018) 

17.2 

22.9 

12  (210) 

24.6 

35.2 

5  (024) 

20.3 

31.8 

13  (240) 

52.1 

67.6 

6  (030) 

34.3 

40.1 

14  (270) 

59.4 

78.8 

7  (060) 

52.8 

71.3 

15  (300) 

71.9 

74.4 

8  (090) 

63.5 

82.8 

16  (330) 

42.4 

49.8 

Mode  Shape  for  Length  L/32,  Area  Ratio  0.7,  and  f=426.2 
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Experiment  19. 


Constriction 


length;  U  32  ,  area  ratio:  0.7  ,  angle  from  driver:45‘ 


439.3  Hz 


1  (000) 


2  (006) 


3  (012) 


4  (018) 


5  (024) 


6  (030) 


7  (060) 


8  (090) 


Vrms 

Corrected 

(mV) 

Data 

43.4 

43.4 

42.0 

43.4 

34.3 

43.8 

31.7 

42.3 

26.7 

41.8 

34.8 

40.7 

20.2 

27.3 

7.0 

9.1 

Mode:  high 


Vrms 

(mV) 


14.6 


24.0 


51.6 


27.9 


19.8 


5.8 


18.8 


30.6 


9  (120) 


10  (150) 


11  (180) 


12  (210) 


13  (240) 


14  (270) 


15  (300) 


16  (330) 


Corrected 

Data 


18.4 


35.8 


43.6 


40.0 


25.7 


Mode  Shape  for  Length  L/32,  Area  Ratio  0.7,  and  f=439.3 


— Theory 
*  Measurement 
Edge  of  Const. 
Driver  Position 


200  250 

Theta  (degrees) 


300  350 


5 


Experiment  20. 


1  Constriction 

732  ,  area  ratio:  0.7  ,  angle  from  driver;45° 

Resonance  frq.:  426.9  Hz 

Q:  125 

Mode:  low  Te 

Corrected 

Data 

^Dcrophone 
#  (0) 

Vrms 

(mV) 

Corrected 

Data 

1  (000) 

8.7 

8.7 

9  (120) 

42.3 

53.1 

2  (006) 

6.5 

6.7 

10  (150) 

23.2 

34.5 

3  (012) 

9.5 

12.1 

11  (180) 

9.6 

8.1 

4  (018) 

13.7 

18.2 

12  (210) 

19.1 

27.4 

5  (024) 

15.7 

24.5 

13  (240) 

38.7 

50.2 

6  (030) 

26.1 

30.5 

14  (270) 

43.6 

57.8 

7  (060) 

39.2 

52.9 

15  (300) 

52.5 

54.3 

8  (090) 

46.7 

60.9 

16  (330) 

30.9 

36.3 

Mode  Shape  for  Length  L/32,  Area  Ratio  0.7,  and  f=426.9 
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